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The Landau-Lifshitz (LL) equation, originally proposed at the macrospin level, is increasingly used 
in Atomistic Spin Dynamic (ASD) models. The models are based on a spin Hamiltonian featuring 
atomic spins of fixed length, with the exchange introduced using the Heisenberg formalism. ASD 
models are proving a powerful approach to the fundamental understanding of ultrafast magnetisation 
dynamics, including the prediction of the thermally induced magnetisation switching phenomenon 
in which the magnetisation is reversed using an ultrafast laser pulse in the absence of an externally 
applied field. The paper outlines the ASD model approach and considers the role and limitations 
of the LL equation in this context. 


I. INTRODUCTION 

Atomistic spin models have a long history, going back 
to the pioneering work of Binder^ and co-workers in the 
1970s. Typically these studies concentrated on the static 
properties of spin systems, particularly using Monte- 
Carlo methods to investigate the order/disorder phase 
transition and finite size effects in magnetic nanoparti¬ 
cles. The atomistic approach proved a powerful tool in 
the study of purely thermodynamic aspects of magnetic 
spin systems. 

However, the use of atomistic models for simulating 
magnetisation dynamics was, until recently, rather lim¬ 
ited and generally based on Monte-Carlo simulations of 
escape over energy barriers. Such an approach is ham¬ 
pered by the fact that the timesteps are generally not 
quantified. Time quantification was attempted by Nowak 
et al.l^, but this is successful only in the strong damping 
regime where the precession can be neglected. The study 
of dynamic phenomena however was intrinsically limited 
until the development of dynamicpEl and stochastic atom¬ 
istic spin models 

In general the model of choice for magnetisation rever¬ 
sal studies is micromagnetics. The history of micromag¬ 
netics starts with a 1935 paper of Landau and Lifshitz 
on the structure of a wall between two antiparallel do¬ 
mains, and several papers by Brown around 1940. A 
detailed treatment of micromagnetism is given by Brown 
in his 1963 boolP^. For many years micromagnetics was 
limited to the use of standard energy minimization ap¬ 
proaches to determine domain structures and classical 
nucleation theory to investigate magnetization reversal 
mechanisms in systems with ideal geometry. Arguably, 
the current interest in micromagnetics arises from the 
availability, from about the mid-1980s onward, of large- 
scale computing power which enabled the study of more 
realistic problems which were more amenable to compar¬ 
ison with experimental data. One important realization 
during this period was the fact that, although micromag¬ 
netics can predict the nucleation fields for the magnetic 


system, due to the coexistence of different energy min¬ 
ima, multiple magnetization reversal paths are possible. 
Thus micromagnetics does not predict necessarily cor¬ 
rectly the state of the system after magnetization rever¬ 
sal. Consequently, a lot of work has gone into the de¬ 
velopment of dynamic approaches which use simulations 
based on the Landau-Lifshitz equation of motion. This 
is probably the technique in most common use today. 


Dynamic calculations using micromagnetics have be¬ 
come ubiquitous, finding applications in fundamental in¬ 
vestigations of reversal dynamics of magnetic materials. 
In addition micromagnetic models are vital to many in¬ 
dustries, including the development of new generations 
of magnetic recording heads and media and permanent 
magnets. However, limitations of the micromagnetic ap¬ 
proach are becoming increasingly apparent. Firstly, as 
magnetic materials become increasingly structured at the 
nanoscale to investigate new physical phenomena and 
create new functionality, the continuum nature of mi¬ 
cromagnetics reaches the limits of validity. Secondly, 
although thermal activation can be introduced into the 
micromagnetic formalism, its applicability is strictly lim¬ 
ited to low temperatures; it is known that micromag netic 
models greatly over-estimate the Curie temperature.l^E^ 
This is a serious limitation in the investigation of ul¬ 
trafast magnetisation dynamics, where temperatures up 
to and beyond can be achieved on the picosecond 
timescale. The natural evolution is toward dynamical 
approaches with atomistic resolution; such models are 
gaining increasing traction in dealing with the physics 
of ultrafast magnetisation processes and practical prob¬ 
lems such as providing an und erstan ding of Heat Assisted 
Magnetic Recording (HAMR).IMJ12I Remarkably, the key 
to Atomistic Spin Dynamic (ASD) models is the use of 
the LL equation at the atomic level. Here we outline the 
basis of ASD models, review some recent simulations of 
ultrafast spin dynamics and consider the physical justi¬ 
fication for the use of the LL equation at the atomistic 
level. 
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II. ATOMISTIC SPIN DYNAMIC MODELS 

The physical basis of the atomistic spin model is the 
localization of unpaired electrons to atomic sites, leading 
to an effective local atomistic magnetic moment, which is 
treated as a classical spin of fixed length. Ab initio calcu¬ 
lations of the electron densitjff^ show that in reality, even 
in itinerant ferromagnets, the spin polarization is well lo¬ 
calized to the atomic sites. Essentially this suggests that 
the bonding electrons are unpolarised, and after taking 
into account the bonding charge the remaining d elec¬ 
trons form a well-defined effective localized moment on 
the atomic sites. Nonetheless the assumption of classical 
spins leads to a fundamental discrepancy with experi¬ 
ments which will be discussed later. 

The basis of ASD models, reviewed by Evans et al.lSl 
is a classical spin Hamilonian based on the Heisenberg 
exchange formalism. The spin Hamiltonian % typically 
has the form: 


'H — Hexc + 'Hani + 'Happ, (1) 

with the terms on the RHS representing respectively 
the exchange, anisotropy and Zeeman terms. The ex¬ 
change term is usually isotropic in spin space and the 
anisotropy term includes energies which are angular de¬ 
pendent. These can arise from crystalline anisotropies 
or magnetostriction and strains. The exchange term can 
also be anisotropic in some situations as will be discussed 
later. 

The exchange energy for a system of interacting atomic 
moments is given by the expression 

^exc = -^JyS,.S, (2) 


where is the exchange interaction between atomic 
sites i and j. Si is a unit vector denoting the local spin 
moment direction and Sj is the spin moment direction 
of neighboring atoms. The unit vectors are taken from 
the actual atomic moment /Xg and given by Si = /Xs/|/Xs|. 
Due to the strong distance dependence of the exchange 
interaction the sum in Eq. © is often truncated to in¬ 
clude nearest neighbors only. This significantly reduces 
the computational effort while being a good approxima¬ 
tion for many materials of interest. In reality the ex- 
cha nge i nteraction can extend to several atomic spac- 
ing ji5 | i6 j, representing hundreds of pairwise interactions. 

In the simplest case the exchange interaction Jij is 
isotropic, meaning that the exchange energy of two spins 
depends only on their relative orientation. In more com¬ 
plex materials, the exchange interaction forms a tensor 
with components: 


JJy — 


Jxx Jxy Jxz 
Jyx Jyy Jyz 

Jzx Jzy Jzz 


( 3 ) 


which is capable of describing anisotropic exchange 
interactions, such as two-ion anisotropjEH and the 


Dzyaloshinskii-Moriya interaction (off-diagonal compo¬ 
nents of the exchange tensor). In the case of tensorial 
exchange, the exchange energy is given by the product: 

= (4) 

We now proceed to consider two important factors in 
the use of ASD, firstly the process of determining, from 
first principles, the parameters of the spin Hamiltonian 
and secondly the introduction of spin dynamics, and the 
implications of the use of the LL equation at the atomistic 
level. 


A. Ab-initio calculation of spin Hamiltonian 
parameters: Multiscale approaches 

The material parameters central to the spin Hamilto¬ 
nian are the exchange interactions, anisotropy energies, 
and the magnitude of the spin moment. These can gen¬ 
erally be found using two routes: (i) experimental mea¬ 
surements, either in a mean-field sense from macroscopic 
quantities such as the Curie point, or microscopically us¬ 
ing neutron scattering, (ii) with a multiscale approach 
using ab initio Density Functional Theory (DFT) calcu¬ 
lations to parameterise the spin model. The ab initio 
approach is often preferable as it removes extrinsic fac¬ 
tors from the parameters such as non-uniformity of an 
experimental sample and also provides a resolution (for 
example many exchange neighbours) which is hard to ob¬ 
tain experimentally even with Neutron scattering. Diffi¬ 
culties can arises such as in magnetic materials involving 
Rare-Earths, where the treatment of the 4/ electrons is 
problematic in DFT. In this case an experimental pa¬ 
rameterization becomes the most practical route, as was 
done by Ostler et al.^for amorphous GdFeCo alloys and 
Evans et al.l^ for Nd 2 Fei 4 B alloys. 

Contemporary ab initio methods enable the calculation 
of a wide range of material properties including ground 
state magnetic properties. So-called ‘beyond DFT’ meth¬ 
ods allow the calculation of even small energy differences, 
providing access to the magnetic crystalline anisotropy 
constants. Standard software packages such as VASP^ 
and SIESTAI^ make such calculations accessible to in¬ 
terested researchers. The calculation of pair wise ex¬ 
change interactions in DFT is somewhat complicated by 
the delocalised nature of the basis sets employed. Pa¬ 
rameterizing a spin Hamiltonian therefore requires map¬ 
ping many different spin configurations onto the atomic 
Hamiltonian. An alternative is to u se sc attering meth¬ 
ods such as Korrin ga-Kohn-RostokeJ^l^^ or linear muf¬ 
fin tin orbital^^Hm These methods are built around the 
atomic sphere approximation which gives a natural map¬ 
ping onto the localised Heisenberg formalism in conjunc¬ 
tion with the magnetic force theoreirP^. 

Connecting the ab-initio and atomistic length scale is 
an important link in the multiscale modelling chain al¬ 
lowing one to include both dynamics and temperature, 
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which has been demonstrated for FePt by Kazantseva 
et alP^. In this section we consider the calculation of 
the temperature dependence of static materials proper¬ 
ties. Mryasov et al.l^ carried out ab-initio calculations 
of exchange and anisotropy of the LiO phase of FePt; an 
important candidate for HAMR media. The aim was to 
investigate the exchange and anisotropy values of bulk 
FePt and to map them onto a classical spin model in 
order to investigate static and dynamic properties. This 
process was complicated by the induced Pt moments aris¬ 
ing from the Fe exchange field. Such an effect is beyond 
the Heisenberg moment of fixed moments, but based on 
ab-initio calculations an effective spin Hamiltonian was 
developed which is dependent only on the Fe degrees of 
freedom. 


^FeP‘ = -^(4S,.S, 






( 5 ) 

The exchange parameters include the effect of Fe-Pt-Fe 
interactions which contributes to both the isotropic ex¬ 
change Jij as well as introducing a two-ion anisotropy, 

Dij , because of the layered ordering of this intermetal- 
lic. This two-ion term is considerably larger than the 
single ion anisotropy ■ The exchange interactions 
Jij are significant over large distances, making numerical 
calculations rather time consuming. Physically it leads to 
strong finite size effectPi^ and in particular leads to devi¬ 
ations of the finite size scaling exponents from the expec¬ 
tations of the nearest neighbour Heisenberg fornP2H2il_ 

In ref. 1271 the thermodynamic properties of FePt were 
investigated using an atomistic model based on the spin 
Hamiltonian given in Eq. ([^. It was shown that the two- 
ion term gives rise to a thermal anisotropy scaling expo¬ 
nent n = 2 { with K{T) cx M‘^{T)) consistent with the 
theory of Callen and CallerP^ and in contrast to the sin¬ 
gle ion anisotropy for which the exponent is n = 3. The 
importance of the atomistic approach is that it is able 
to calculate the exact exponent arising from the specific 
material parameters of FePt. This resulted in n = 2.1, 
the non-integer value reflecting the relative importance of 
the single- and two- site anisotropy terms; a value i n good 
agreement with experiments on FePt nanoparticleaSSM! 


B. Langevin Dynamics and the LL equation at the 
atomic level 


bath, with the Landau-Lifshitz equation augmented by 
white-noise fields, effectively producing the Stochastic 
(Langevin) equation of the problem. The approach is 
to determine local fluctuating fields using the fluctuation 
dissipation theorem and to require the equilibrium distri¬ 
bution function of the orientations of the magnetization 
to coincide with the Boltzmann distribution. 

In order to yield the Boltzmann equilibrium distribu¬ 
tion, the stochastic LL equation should be interpreted as 
a Stratonovich vector stochastic differential equatioiP^. 
This is integrated by a suitable choice of the numerical in¬ 
tegration scheme, most usually that of Heun. Care must 
be taken that the spin moments remain of unit length 
and for a non-conservative scheme such as Heun, an ex¬ 
plicit renormalization of the length at each timestep is 
required for the Stratonovich solutiorPH. The integration 
of the stochastic LL equation is discussed in detail in 
ref. UH 

Consequently, the basis of ASD for a set of coupled 
spins is the integration of the stochastic Landau-Lifshitz 
equation for each localized magnetic moment S^: 

S, = 7[S, X H,] - 7a[S, x [S, x H,]] (6) 

Here = ^j(t) — dH/dSi is the local effective field 
which includes Zeeman, exchange, anisotropy and mag¬ 
netostatic contribution, augmented with by a stochas¬ 
tic term (which appears like a field). It is defined 
through the correlators: 

m)) = 0 , ( 7 ) 

Here T is the temperature of the heat bath, 7 is the gy- 
romagnetic ratio, /ig is the magnetic moment, a is the 
parameter describing the coupling strength to the heat 
bath, 77 and v are Cartesian components. The basis of 
this equation is the separation of timescales, assuming 
that the bath (phonon or electron system) is much faster 
than the spin system. Consequently, the fluctuation- 
dissipation theorem can be applied to derive the equi- 


we consider the applicability of the Langevin Dynamic 
approach using the LL equation. However, prior to this 
we give some examples of the success of ASD models in 
developing an understanding of the thermodynamic as¬ 
pects of ultrafast magnetisation processes, including the 
prediction of a novel ‘linear’ reversal mechanism. 


librium white noise properties of Eq. Q. In section IV 


The important step forward in the use of atomistic 
models is the introduction of Langevin Dynamics, al¬ 
lowing modelling of the dynamical response of the mag¬ 
netisation to temperature change^^. The approach is 
based on the introduction of thermal fluctuations for a 
single particle developed by BrowrP^. The theoretical 
basis is the classical theory of Brownian motion which 
accounts for the departure from thermal equilibrium due 
to the energy interchange between a particle and its heat 


III. ATOMISTIC MODELS OF ULTRAFAST 
SPIN DYNAMICS 

Advanced models are required in order to understand 
ultrafast magnetisation processes. Interest in this area 
has developed rapidly since the experimental demonstra- 
tiorP^ that the magnetisation of Ni can be reduced by 
laser heating on a timescale of Ips. Experimentally the 
measurements are made using a pump-probe process. A 
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high energy femtosecond laser is used to heat the mate¬ 
rial (pump), and the magnetic response is measured us¬ 
ing MOKE with a low energy probe beam split off from 
the pump. This experiment gives time resolved measure¬ 
ments of the magnetic response following pulsed laser 
heating. Such experiments are extremely challenging in 
terms of understanding the physics of ultrafast magneti¬ 
sation processes and damping mechanisms. 

A more recent development was the experimental 
demonstration by Stanciu et al.® of optically induced 
magnetisation reversal in the amorphous ferrimagnet 
GdFeCo. Using circularly polarised ultrafast laser pulses, 
Stanciu et al. showed magnetisation reversal to be de¬ 
pendent on the chirality of the laser pulse. This was 
interpreted as arising from a large, laser-generated, opto- 
magnetic field (estimated as large as 20T) possibly origi¬ 
nating from the Inverse Faraday Effect. The reversal was 
explained^ using atomistic and Landau-Lifshitz- Bloch 
(LLB) macrospin simulations (for a review of the LLB 
equation see^ as arising from this large optomagnetic 
field together with purely thermodynamic contribution 
which initiates switching via the so-called linear reversal 
mechanisirP^. Linear reversal is an important prediction 
of the atomistic model and property of the LLB equa¬ 
tion. Essentially it involves a collapse of the magneti¬ 
sation to zero and subsequently a switched polarisation 
in a reversing field. Linear reversal sets in at a critical 
temperature T* related to the ratio of the longitudinal 
and transverse susceptibilities^. Importantly, this leads 
to ultrafast magnetisation reversal since the timescale is 
governed by the longitudinal relaxation time which is in 
the order of hundreds of femtoseconds. 

Originally optomagnetic switching was only observed 
in various ferrimagnetic structures but it can, in prin¬ 
ciple, occur in ferromagnets as well. The experimental 
work by Lambert et al)^ found that thin films of ferro¬ 
magnetic Co/Pt layers and granular FePt showed mag¬ 
netisation switching after repeated excitation by a cir¬ 
cularly polarised femtosecond laser. Fig. shows an 
atomistic model simulation of the response of the mag¬ 
netisation of FePt to the field and temperature pulses 
associated with a femtosecond laser pulse. The calcula¬ 
tions used the spin Hamiltonian of Mryasov et given 
in Eq. ([^ with a optomagnetic field of SOT that lasts for 
1 ps after laser excitation. The dynamic response is cal¬ 
culated using Langevin Dynamics driven by the electron 
temperature evolved using a 2-temperature modeP^. The 
reversal depends critically on the laser power, as shown in 
Fig. 0 At low laser power. Fig. [^a), demagnetisation is 
not complete and the reversal proceeds via the usual pre- 
cessional route, albeit over an energy barrier reduced due 
to the quenching of the anisotropy. Complete reversal 
proceeds over many picoseconds. At elevated tempera¬ 
ture, Fig. Bt). switching of the total magnetisation pro¬ 
ceeds by a process involving no macroscopic precession; 
the linear reversal mode. Importantly, switching via the 
linear reversal mode occurs on a timescale of the longi¬ 
tudinal relaxation time of the magnetisation (^hundreds 
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FIG. 1. The temporal evolution of the magnetization in FePt 
after a femtosecond laser pulse simulated using the atomistic 
model described in section II. The model uses the Hamiltonian 
expressed in Eq. [^with a damping, a = 0.1, and fis = 3.23/rB. 
The system starts at room temperature before excitation by 
a 100 fs laser pulse incident at t=0 ps with fiuences of (a) 
6 mJ/cm^ (b) 12 mj/cm^ and (c) 18 mj/cm^. A reversing 
optomagnetic field also occurs with the laser pulse; it has a 
square pulse shape from t=0 ps to t = 1 ps with a magnitude 
of 30 T. 

of fs). At higher laser power, Fig. [^c), reversal occurs 
but the magnetisation is destroyed by the elevated tem¬ 
perature. Vahaplar et al.l^ show that optically induced 
switching proceeds above a critical temperature which is 
sufficient to excite linear reversal but not so large as to 
demagnetise the system. 

It was shown by Vahaplar et al.l^ that extremely high 
optically-induced fields (tens of Tesla) needed to be in¬ 
voked to trigger the all-optical reversal. The possible 
origin of such fields remains a matter of debate. How¬ 
ever, it can be shown that the fields may arise from inter¬ 
atomic exchange forces. This interpretation begins with 
the observation by Radu et al.l^ that the RE and TM 
sublattices demagnetise at different rates, even though 
strongly coupled through intersublattice exchange forces. 
The measurements were made using XMCD to provide 
the element specific magnetisation dynamics. The exper¬ 
imental observations were supported by atomistic model 




























5 


calculations, which verified both the differential sublat¬ 
tice dynamics and the existence of an intriguing transient 
ferromagnetic-like state (TFMLS) which is created by the 
reversal of the TM spins into the RE spin direction. The 
TFMLS exists for around SOOfs and is associated with 
the magnetisation reversal process. Further investigation 
led to the astonishing prediction by Ostler et al.l^, us¬ 
ing atomistic model calculations that switching occurred 
in the absence of any external symmetry breaking field. 
In ref. im this remarkable prediction was verified experi¬ 
mentally. Interestingly, thermally induced magnetisation 
switching (TIMS) allows a re-interpretation of optically- 
induced magnetisation reversal. Rather than invoking 
large fields of opto-magnetic origin, Khorsand et. al.SSl 
gave experimental evidence in favour of TIMS as the 
main switching mechanism in GdFeCo and attribute the 
helicity dependence of the laser excitation on the dichroic 
effect, i.e. the dependence of the absorption of energy 
from the laser pulse on the chirality of the laser light. 

The phenomenon of TIMS has been furt her in ves- 
tigated theoretically using macrospin model^^^*^ and 
atomistic approache^^. Importantly, refs. [SUI and IM1 
show that TIMS arises from angular momentum transfer 
between the RE and TM sublattices mediated by the es¬ 
tablishment of a two-magnon bound state. The study in 
ref. [FT] involved detailed calculations of the magnon band 
structure, which has two branches with properties which 
strongly depend on the material composition. Calcula¬ 
tions showed that transfer of energy between the modes, 
resulting from laser excitation, was the physical origin of 
the TIMS phenomenon. Further, the calculations demon¬ 
strated a window for TIMS within a certain range of al¬ 
loy concentrations. Specifically, for low RE concentra¬ 
tion essentially a uniform FM mode is excited. With 
increasing RE concentration, the optical mode becomes 
accessible, leading to TIMS. At higher RE concentrations 
there develops a large band gap which precludes the an¬ 
gular momentum transfer between sublattices, at which 
point TIMS cannot be excited. This prediction is in good 
agreement with experiment. 

The results of ref. [5T] allow the definition of design rules 
for TIMS. Specifically, 

1. The existence of two sublattices with differential 
demagnetisation dynamics 

2 . Antiferromagnetic coupling between the sublattices 

3. Reversal of the dominant magnetic sublattice which 
stabilises the switched magnetisation direction dur¬ 
ing the cooling phase following the laser pulse. 

This gives rise to the expectation that synthetic ferri- 
magnetic (SFIM) structures consisting of two ferromag¬ 
netic layers separated by an element such as Ru or Ir, to 
establish AF coupling between the layers, would exhibit 
TIMS. This has been demonstrated numerically by Evans 
et al.^in bilayers of Fe and FePt coupled by an exchange 
separation layer assumed to promote AF coupling be¬ 
tween the FM layers. Importantly, this extends the TIMS 


phenomenon beyond its initial prediction and discovery 
in amorphous ferrrimagnetd^ to designed materials with 
high anisotropy and the avoidance of RE components. 
This is of importance both in terms of fundamental un¬ 
derstanding of ultrafast magnetisation processes and of 
applications in information storage. In the latter context 
we note that Evans et al.l^ demonstrated the requirement 
for large write fields in the magnetic recording process; 
not simply to ensure magnetisation reversal, but also that 
there is no back-switching of the magnetisation, which 
would lead to a limiting source of noise. Ostler et al.l33 
show that extremely large fields ^ tens of Tesla are nec¬ 
essary to oppose the formation of the TFMLS. This is 
consistent with the physical origin of TIMS in the exci¬ 
tation of the two-magnon bound state, which naturally 
introduces fields arising from the exchange interaction. 

Clearly, TIMS is an important prediction using a 
model based on LL dynamics at the atomistic level and 
shows the LL equation to be a remarkable piece of insight 
which finds application at time- and length- scales signif¬ 
icantly beyond its original conception. In the following 
we consider the underlying physics of the LL equation at 
the atomistic level, distinguishing between the origin of 
the precession and damping terms. 

IV. ANALYSIS OF THE ATOMISTIC 
LANGEVIN DYNAMIC APPROACH USING THE 
LL EQUATION 

The precession term in the LL equation derives directly 
from the equation of motion 

ihj^{M) = {[M,H]), (8) 

with H = —M.B{t), which leads directly to the lar- 
mor precession term. The damping term is more diSi- 
cult to justify at the atomistic level. Pragmatically one 
can make the case that some term coupling the spin to 
the heat bath is necessary to ensure eventual thermal 
equilibrium with the heat bath. In this spirit, the LL 
damping term is the simplest form capable of including 
this key physical requirement. Some work has been car¬ 
ried out to investigate spin dynamics at the quantum 
level. In particular, Cappelle and Cyorffji^ have inves¬ 
tigated magnetisation dynamics using time dependent 
density functional theory. They construct a gradient- 
dependent density functional which is then used to derive 
the phenomenological LL form of damping directly from 
first principles. The LL equation and its damping term 
can also be derived from th e q uantum mechanics using 
the density matrix formalisni^, an approach recently re¬ 
visited by Weiseil^. The coupling of the spin to the heat 
bath is expected to be material dependent; essentially 
it represents the relevant energy transfer channel which 
might arise from interaction between the spin and con¬ 
duction electron or spin/lattice interactions and should 
be studied at the quantum level. 
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The introduction of thermal fluctuations into the 
atomistic model formalism is achieved via the Langevin 
Dynamic approach. As mentioned previously, the as¬ 
sumption is made that the instantaneous random fields, 
by which the thermal fluctuations are introduced, are 
uncorrelated in time and space, giving rise to white 
noise. The basis of this formalism is the separation of 
timescales, assuming that the bath (phonon or electron 
system) is much faster than the spin system. In this case, 
the bath degrees of freedom can be averaged out and re¬ 
placed by a stochastic field with white noise correlation 
functions. In sense this is perhaps the most questionable 
part of the application of the atomistic model in ultrafast 
magnetism, which involves phenomena on the timescale 
of tens of femtoseconds, where spatial and temporal noise 
correlations might be expected. 

Investigation of this phenomenon is the province of 
models simulating the thermal behaviour of the heat bath 
and coupling this to the spin system. In this case only 
the precession term of the LL equation is retained, with 
the damping included via a term based on the underly¬ 
ing mechanism coupling the spin to the heat bath. This 
was treated in a generic way by Atxitia et al.^ using an 
approach developed by Miyazaki and SekP^ and gener¬ 
alising this to multi-spin systems. The noise takes the 
form of an Ornstein-Uhlenbeck proces^^ introducing a 
characteristic correlation time Tc- The spin is coupled 
locally to the bath which is connected to a thermostat as 
follows 

Si = 7[Si X (Hi -h T]i)], 

m = -—{vt- xSj) + Ri (9) 

Tc 

with the fluctuation-dissipation theorem for the 
bath variable: (Ri(i)) = 0 ; (Ri(t)Rj(t')) = 

{2xkBT/Tc)Sij6{t — t'). The parameter x describes 
the coupling of the bath variable to the spin. The 
precession term in the first equation of the set ^ has 
the same form as in the Eq. §. However, the damping 
is now described by the second equation in this set. 
In the second equation in Eqs. the bath variable 
adjusts to the direction of the spin due to the interaction 
with it. In the limit Tc —>■ 0 the stochastic LL equation 
is recovered. This also provides a relation between 
the damping and the coupling constants as a = qxT;, 
giving a more precise physical sense to the LL damping 
constant at the atomistic level. 

In this approach, the phenomenological LL damping 
parameter is substituted by two unknown parameters: 
the correlation time Tc and the coupling constant x- Sev¬ 
eral processes may be important in determining these 
constants, for example, the spin-orbit coupling, momen¬ 
tum relaxation, scattering rate and dephasing time of 
conduction electrons. As in the LL approach, these pa¬ 
rameters will be material specific and their physical ori¬ 
gins should be clarified on the basis of first-principle ap¬ 
proaches. The effect of the correlation time on the ultra¬ 
fast demagnetisation process was investigated in ref. |S3 


For correlation time Tc < 1 fs the correlated approach 
gives the same results as the standard Langevin Dynam¬ 
ics with white noise However, in the range Tc — 10 — lOOfs 
the correlations were found to give a dramatic increase 
of the longitudinal relaxation time. The effect is less pro¬ 
nounced at higher temperature since in this case the tem¬ 
perature contributes to the loss of correlations. Calcu¬ 
lations based on the Langevin Dynamic approach gener¬ 
ally give reasonable values for the longitudinal relaxation 
time in comparison with experiment, which suggests that 
experimental correlation times are on the order of 10 fs 
or less; greater values would have an appreciable effect 
on the observed rates of demagnetisation. 

The interactions of the spins with the lattice system 
also provides an energy channel for the fluctuation and 
dissipation. In conventional atomistic models the lattice 
is fixed and so transfer to and from the lattice is han¬ 
dled phenomenologically by the fluctuation and damp¬ 
ing terms in the stochastic LL equation. By introducing 
the motion of the atoms to the model the energy trans¬ 
fer can be modelled directly without the need for the 
phenomenological terms. Recently models such as this 
have been developed to investigate a variety of systems 
where the spin-lattice effects are important. Ma et 
have extensively developed a spin-lattice model of Fe that 
utilises the dependence of the exchange on the atomic 
separation as the coupling between the spins and lat¬ 
tice. However in this case both the spins and lattice use 
a Langevin thermostat to maintain a constant tempera¬ 
ture. Through this model both systems act as a thermal 
reservoir from which instantaneous spin and lattice tem¬ 
peratures can be extracted which can then be dynami¬ 
cally linked to the electron temperature thus representing 
a dynamic three temperature model.^ Using this model, 
termed Spin-Lattice-Electron Dynamics (SLED) by Ma 
et al. , experimental ultrafast laser induced magnetisa¬ 
tion dynamics can be reproduced. 

Whilst the exchange interaction can couple the spins 
and lattice further coupling terms have also been inves¬ 
tigated. Karakurt et implemented a spin-lattice 

model where the exchange is constant but introduce a 
specific coupling term of the form: 

Uc = -cs, ■ r,, (10) 

Where is the separation of the atoms and 

C is a parameter to control the strength of the coupling. 
Using this Karakurt et al. demonstrated that this cou¬ 
pling causes a damping of the uniform precession mode. 
Beaujouan et propose a different type of coupling 
based on a two site anisotropy which takes the form of 
pseudo-dipole interaction. 

He = -K{r,j) ^(fij • Si) (r„- • S^-) - (H) 

In this case the coupling strength, K, depends on the sep¬ 
aration of the atoms and requires specific parametrisation 
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from ab initio. As discussed in the literature this form of 
coupling arises from the spin-orbit interaction of the elec¬ 
trons and thus is more physically justifiable but is still 
not exact. With this Beaujouan et al. are able to show 
that energy can be transferred between the systems and 
an equilibrium temperature is obtained. It is clear that 
by incorporating lattice dynamics into the ASD model 
various effects which are treated phenomenologically are 
present. However these still require a high level of em¬ 
pirical parametrisation for both the atomic bonding and 
the spin-lattice coupling. 

Finally, we consider briefly one further aspect of the 
use of fixed spin models such as the LL equation, specifi¬ 
cally the classical nature of the spin. This leads to a dis¬ 
parity between the simulated and experimental tempera¬ 
ture dependent magnetization curve^^. At the macro¬ 
scopic level the temperature dependent magnetization 
is well fitted by the phenomenological equation pro¬ 
posed by Kuz’mirP^. However, the Kuz’min equation 
merely describes the form of the curve with little re¬ 
lation to the microscopic interactions within the mate¬ 
rial which determine fundamental properties such as the 
Curie temperature. Ideally one would perform 3D quan¬ 
tum Monte Carlo simulations^, but, although this is 
possible for small numbers of atoms, for larger ensem¬ 
bles the multiscale approach using atomistic models pa¬ 
rameterized with ab-initio information remains the only 
feasible approach to connect the quantum and thermo¬ 
dynamic worlds. Evans et al.^^ have proposed a scaling 
approach which maps classical to quantum spin models. 
The scaling recognises that, although the classical treat¬ 
ment finds the correct magnon energies, the distinction 
between classical and quantum models results from the 
particular statistical properties of each approach. While 
quantum systems obey Bose-Einstein statistics, leading 
to the Bloch law at low temperatures, the classi¬ 
cal Boltzmann statistics gives rise to a finite slope of the 
magnetisation at low temperatures. In ref. the ex¬ 
istence of a simple relation between classical and quan¬ 
tum temperature-dependent magnetization at low tem¬ 
peratures is demonstrated. The temperature-dependent 
magnetization is represented in the simplest form arising 
from a straightforward interpolation of the Bloch law and 
critical behaviour given by the Curie-Bloch equation 

m(T) = (l-r“)^ (12) 


where a is an empirical constant and /3 « 1/3 is the crit¬ 
ical exponent. Evans et al.^ then use the classical spin 
model simulations to determine the critical exponent /3 
and then find a by fitting the classical model predictions 
to experimental data. This leads to a mapping from a 
‘simulation temperature’ to the real temperature. It was 
shown that this approach gives excellent agreement with 
experimenlP^ for the demagnetisation of Ni following an 
ultrafast laser pulse. 

V. CONCLUSION 


The use of the LL equation in ASD models of mag¬ 
netic materials has been described. The introduction of 
the LL equation, in its stochastic form, is the basis of a 
powerful approach to ultrafast sin dynamics. In partic¬ 
ular, ASD simulations demonstrate the important ther¬ 
modynamic contribution to laser-induced ultrafast pro¬ 
cesses. The models demonstrate a new, so-called linear 
reversal mechanism which is the path to ultrafast rever¬ 
sal. Also predicted is the phenomenon of TIMS, which 
is currently under extensive investigation and holds the 
promise of application in future devices requiring fast 
switching. The success of the LL equation in this frame¬ 
work is remarkable. While further investigations of the 
energy transfer mechanisms at the quantum level should 
be carried out to improve the physical understanding 
of damping, ASD methods based on the LL equations 
are likely to have an important role in understanding 
the physics of magnetic phenomena, not only at short 
timescales and elevated temperatures, but also on length- 
scales where the micromagnetic formalism is not appro¬ 
priate. 
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